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Abstract 

Certain spin Hamiltonians that give rise to tunnel splittings that are viewed 
in terms of interfering instanton trajectories, are restudied using a discrete 
WKB method, that is more elementary, and also yields wavefunctions and 
preexponential factors for the splittings. A novel turning point inside the 
classically forbidden region is analysed, and a general formula is obtained for 
the splittings. The result is appled to the Fes system. A previous result for 
the oscillation of the ground state splitting with external magnetic field is 
extended to higher levels. 
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The magnetic properties of the molecular cluster [(tacn) 6 Fe802(OH)i 2 ] 8+ (or just Feg for 
short) are governed by a Hamiltonian |l 

H = -k 2 J 2 z + {k 1 -k 2 )J*-g f i B J-H, (1) 

where J is a dimensionless spin operator, H is an externally applied magnetic field, J = 10, 
k\ ~ 0.33 K, and k 2 ~ 0.22 K. The zero-field Hamiltonian has biaxial symmetry with easy, 
medium, and hard axes along z, y, and x respectively J2|. Very recently, Wolfgang and 
Sessoli P] have seen a new effect in this system, viz., an oscillation in the Landau-Zener 
transition rate between Zeeman levels, as a function of the applied field along x. These 
oscillations reflect an oscillation in the underlying tunneling matrix element between the 
levels in question, and are in this author's view, the only unambiguous evidence to date for 
quantum tunneling of a spin of such a large size in a solid state system. This effect is not 
seen, e.g., in a closely related Mn 12 cluster. 

Oscillations as a function of H x in the ground state tunnel splitting A of the Hamiltonian 
(|I|) were in fact predicted earlier [[|], on the basis of an instanton calculation. For H||x, there 
are two instantons, with a complex action differing by a Berry phase that sweeps through 
odd multiples of 7r as H x is varied, leading to a complete quenching of tunneling. In this 
view the effect arises from destructive intereference between spin trajectories |J. A different 
perspective was provided in Ref. |§, by noting that 7i is invariant under a 180° rotation 
about x when H z = 0. The oscillation is due to a symmetry- allowed crossing of levels with 
different parity under the rotation. 

Although the oscillations with H x are easily seen by direct numerical diagonalization of 
H (see Fig. 1 of 0], or Fig. 4 of 0), it is of interest to understand these features analytically. 
Since the spin J is large, it is natural to use the semiclassical, or J — > oo approximation. To 
some extent, the instanton method already does this. In this paper, we will try and make 
further progress using a discrete WKB method pj-|ilfl, which has several advantages. First, 



it is very difficult to find the next-to-leading terms in the J — > oo asymptotic expressions for 
various physical quantities (such as A) using instantons. Second, Wernsdorfer and Sessoli 
also see an oscillatory rate in the presence of a dc field along z that is such as to align the 
ground level in one well with an excited level in the other. The splitting is now never perfectly 
quenched as the symmetry of 7i is destroyed. We do not (although others may) know how 
to solve this problem with instantons. Third, with an eye to the future, the method provides 
wavefunctions in addition to energies, which may be used to calculate matrix elements of 
various perturbations @,|ll|] to Eq. ([!]) that are present in the actual physical system, and 



thus study their influence. 

Specifically, we will derive a general result [see Eq. (|T8|)1 for the tunnel splitting between 
degenerate pairs of levels in a symmetric problem when oscillations are present. Our result 
is expressed in terms of two action integrals. In the course of doing this, we will encounter 
a novel feature that does not arise in previous discrete WKB studies, namely, a turning 
point in the classically forbidden region! We will apply our result to the hamiltonian ([I]) for 
H z = 0, focussing in detail on quenching fields for ground and excited state splittings. The 
results for the latter are new. The study of imperfect quenching of tunneling that occurs 
when H z ^ 0, because the potential is then asymmetric, is much more involved, and will be 
published separately. 

Let us first briefly review the discrete WKB formalism M. The starting point is to write 
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Schrodinger's equation in the J z basis. Let TL\ip) = E\ip), J z \m) = m\m), (m\ip) = C, 
(m\H\m) = w m , and (m\H\rn!) = t mt7n > (m ^ m'). Then we have 



tm,n,C n + W m C m — EC m . (2) 

We assume that the matrix t m „ is real and symmetric, t m ^ n = t n>m . In the present problem, 
we need matrix elements that are off-diagonal by 1 (i m , m ±i) and by 2 (i m , m ±2)- This makes 
Eq. (H|) a recursion relation involving five terms, as opposed to three terms in previous work. 
The physical idea is to view Eq. (TJ) as a tight-binding model for an electron hopping on a 
one-dimensional lattice, and use the approximation of semiclassical electron dynamics. This 
would be exact if the matrix elements of 7i were constant with m, and will be systematically 
justifiable if they are slowly varying with m. Formally, the latter means that we can find 
functions w(m), ti(m), and t2{m) of a continuous variable m, such that on the discrete 
eigenset of J z , 

w(m) = w m , (3) 

t a (rn) (im,m+a tm,m— a)/2, Ot 1, 2, (4) 

and further, that if m/J is regarded as a quantity of order J°, then w(m) = dw/dm = 
0(w(m)/ J), with similar restrictions on ti(m) and t2{m). For Eq. ([!]), these conditions are 
met if J > 1. 

If w m , t m>m ±i, and t m ^ m ±2 were constant, the eigenstates of TC would be states with 
C m = e tqm , and E = w + 2t!Cosg + 2t 2 cos(2g). Now we seek a solution in the form 
C m = e^ m ) with $ = $ + $i + $ 2 + • • •, where $ n = 0(J 1_n ), and $ n = 0($ n /J). 
Then, one can show H] that up to terms of order J° in $, the solution is given as linear 
combinations of the form 

C m ~ exp (i [ q(m')dm'\ (5) 

Jv(m) ^ J ' 

where q(m) is a local wavevector that obeys the eikonal or Hamilton- Jacobi equation, 

E = w(m) + 2ti(m) cosg + 2t 2 (m) cos(2g) = Ti, sc (q, m), (6) 

and u(m) is the associated semiclassical electron velocity, which obeys the transport equation 

v(m) = dH S c/dq = —2sxnq(m)(ti(m) + 4t 2 (?n) cos q(m)). (7) 

To talk of tunneling, we must first understand the classically allowed and forbidden 
regions in the m space. As a function of q for fixed m, the semiclassical Hamiltonian 
?i sc (q,m) can be viewed as a band energy curve, and its minimum and maximum values 
define local band-edge functions U±(m). The classically accessible region for any energy E 
is thus defined by £7_(m) < E < U + (m). [The first consequence of having five terms in the 
recursion relation shows up here. In the three term case, the band edges always occur at 
q = or 7r. Now, they can occur at values other than these if |ti(m)/4t 2 (^)| < 1. These 
functions are sketched in Fig. 1 for Eq. ([I]) with H z = 0. The minimum, [/_, is attained at 
q = for \m\ > m*, and at q ^ for \m\ < m*. The curve U-(m) is smooth at m — ±m*, 
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and the formula for m* is unimportant.] Thus, for the energies E a and Ej, drawn in Fig. 1, 
the central region is classically forbidden and allowed, respectively. We will focus on states 
of the first type in what follows. 

The next step is to derive a generalization of Herring's formula for the tunnel splitting A 
for a pair of levels whose mean energy is E. Proceeding in exact analogy with Ref. JE2J (see 



also ||(c,d) or fllPl), we consider a solution C m to Eq. (g) with energy E, that is (a) localized 
in the left well of U-(m), and decays away from that well everywhere including the region 
near m — 0, and (b) normalized to unit total probability. The behavior of this solution near 
the right well need not be specified or examined too closely. Up to an irrelevant over all 
sign, we find 



A 
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toxoid - C_i) + t 02 C (C 2 - C_ 2 ) + t^iCl - Cl x )] , integer J 



2 t_i i (Ci - C\ +4 t_a i CiCs -C_iC_a , half-integer J. 

2'2\ 5 — 2' 2 ' 2 V 2 2 2 2 / 



(8) 



To apply Eq. @, we must find C m in the central region. In principle the procedure is 
straightforward, and follows conventional WKB. We first find C m in the allowed region, near 
—mo, and then use connection formulas to extend it into the forbidden region. For a three 
term recursion relation, this is done in Ref. [ 10| . In the present case, we encounter a new 



difficulty To see this, we consider points at which v{m) vanishes. At all such points, which 
may be called turning points, the solution @ diverges, indicating a breakdown of the WKB 
approximation. Let us now consider a point strictly inside the classically allowed region in 
the E-m plane. At such a point q is not an extremum of 7i sc for fixed m, i.e., v(m) ^ 0. It 
is a simple corollary that the points E = U±(m) are turning points, corresponding to q — 0, 
7r, or cos _1 (— ti/At 2 ). These turning points are of the same physical character as those in 
conventional WKB, and the q = or it ones are the only ones that arise with a three-term 
recursion relation. 

For our five-term recursion, however, v(m) can also vanish if cosg = —ti/4t 2 , even though 
E 7^ U±{m). To see how this can happen, we solve Eq. (H) to get 



. . -tx(m) ± t?(m) - 4t 2 fm) / (m) 1/2 

COSfllm = ; r . 

yV ; 4t 2 (m) 



(9) 



where f(m) = w(m)—2t2(m)—E. Thus, such a turning point may arise when the discrimiant 
of the quadratic equation for cosg(m) vanishes. Since, by exclusion, such points must 
necessarily lie in a classically forbidden region, where q(m) is not real, it follows that they 
can only arise in problems where \ti(m)/4t 2 (m)\ > 1 for some m. This fact and Eq. (0) 
then imply that at such a point, cosg changes from real to complex, i.e., q changes from 
pure imaginary to complex, and the wavefunction accordingly changes from an exponential 
decay with one sign to a decay with an oscillating sign. 

Since WKB breaks down at the forbidden region turning points, we need connection 
formulas at these points just as for ordinary ones ||. We will publish the derivation of these 
formulas elsewhere, and here we only give the result. Let the discriminant in Eq. (Q) vanish 
at m = m c , and let cosg be real for m < m c , and complex for m > m c . It is convenient 
to define q(m) = iK,(m) with k > in the region m < m c , and to write s(m) = —iv(m) 
everywhere. (This definition renders s(m) > for m > m c .) We consider the decaying 
WKB solution in the region m < m c : 
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C m = — exp | — I K(m')dmf ] , m < m c , (10) 
2\ s(m) 



where A is chosen to be real. For m > m c , we must consider linear combinations of the type 
(H), with two choices for q(m) which we write as 

Qi,2{ m ) — iK(m) ± x{ m )- (H) 

For the solution to continue decaying, we must still have k > 0, and we also choose \ > 0. 
Then, both n(m) and x{ m ) have a kink at m — m c . We further define Si^m) = —iv^qi^im)) 
via Eq. (0), so that s 2 = s*. The WKB solution which connects to (0) is then given by 

A (7 \ 

C m = Re = exp i / qi{m')dm! , m > m c . (12) 

Note that this is explicitly real, as the reality of Eq. (H) requires. Also, Eqs. fllTf ) and (]T2|) 
only hold for \m c — m\ ^> J 1 / 3 . The connection formula for the growing solution is similar, 
but is not needed for our present purpose. 

The result (|8|) is exact, but does not reveal the physically important barrier penetration 
factor. To remedy this, we substitute Eq. ([12]) in Eq. We consider a situation as in 
Fig. 1, with minima in U-(m) at ±mo, and forbidden region turning points at ±mi. The 



key, clearly, is to simplify Eq. (|12" ) in the region \m\ < mi. To this end, we substitute 
Eq. ( |TTD for q in Eq. (^) and separate the real and imaginary parts. This yields 

cosh k cos x = —ti/4:t2, (13) 
sinh k sin X = (4* 2 / - *? ) 1/2 /4t 2 • (14) 

Using these results, it follows that 

Si = 8^2 {m) sinh n{m) sin x( m ) sm <?i( m )- (15) 

We now specialize to the case of integer J; the other is similarly analysed, and yields the 
same result, Eq. flTTD below. For Eq. (|8p, we only need C m for m = 0, ±1, and ±2. The 
variations in k, x, ti, t2, and q± between these points may be ignored as they are of of order 
J~ 2 . Hence, to sufficient accuracy one may write (for \m\ < 2), 

i(n+mq 10 ) 

C m = ReA 2 , (16) 

V sm q w 

where f2 = J_ m qi(m')dm', and A 2 = (8t2o sinh kq sin xo)~ 1 ^ 2 A. The suffix denotes quanti- 
ties evaluated at m = 0; thus gio = qi(0), k = k(0), etc. To the same accuracy as Eq. fll6|) 
one may write t i = ^i(0), and to2 = t-i,\ = ^(0) in Eq. (||). If we use Eq. (|13|), and write 
every thing in terms of t2o, kq, and xo-, then a certain amount of algebra leads to 

A = -A 2 (e 2in + e- 2in ') 

17 \ I 7 \ 

= v4 2 exp — / K,(m')dm' cos / x( m ')dm' . (17) 

\ —mi / \ —mi / 
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The cosine factor clearly shows the possibility of oscillations. 

The next step is to match the WKB wavefunction (10) in the ordinary decaying region 



to the wavefunction in the allowed region. It is plain that A will contain an additional 
barrier penetration factor exp(— f™* K(m)dm), where ±m t are the ordinary turning points. 
Omitting the details of the calculation, which are very much like those in conventional WKB, 
we find that for the nth pair of levels, provided n <C J, 

A = — -g„exp | — I K(m')dm' ] cos | / x( m ')dm' ) , (18) 

where c^o is the small oscillation frequency in the wells near ±m , and g n = \f2iin n e~ n (with 
n — n + |). It need hardly be said that mt and mi depend on the energy and hence on n. 




Equation fll8|) is a general formula for the splitting in the presence of interference effects. 
As opposed to an "exponentially accurate" calculation which gives an asymptotically correct 
result for In A as J — > oo, it is correct for A itself. 

We now apply Eq. fll8|) to the Hamiltonian Eq. (H). The problem is now merely one 
of quadrature, so we will focus only on the cosine factor. (The full expression for the non- 
oscillatory part of A including the exact prefactor is exceedingly lengthy and unilluminating. 
A partial result for the WKB exponent, or Gamow factor may be found in Ref. ||13|| .) In doing 



the quadratures, the first step is to find w(m), t a (m), etc. Here, any function that reproduces 
the first two terms in a series in 1/ J is adequate, since Eq. (^) represents only the two leading 
terms in $(m). We define J = J+h, fi = m/J, H c = 2k 1 J/gfi B , h x = JH X /JH C , A = k 2 /ki, 
and measure all energies in units of k\ J 2 . Then, t\ = —h x (l — At 2 ) 1 / 2 , 4£ 2 = (1 — A)(l — /i 2 ), 
2w = (1 + A)(l — /i 2 ), and cu = 2[A(1 — hi)} 1 / 2 /.!. Further, let us denote the argument of 
the cosine in Eq. (|18D by A. The turning point fii(E) is given by 

o h 2 E 

,^E) = 1 - ^ - -. (19) 

Secondly, from Eq. (|T4"D, we see that x ~ (^i(E) ~ A 4 ) 1 ^ 2 ; so that to relative order 1/J, we 
may write 



M10 / „ 

,9x 
BE 



A = 2j| [x(fi,E = 0)+E 




(20) 



where //io = /ii(0), and E = (n + |)^o f° r the nth pair of levels. At E = 0, we have 
cos 2 x = [(1 — /"io)/(l — A* 2 )]) dx/dE = — cotx/2(l — /i 2 — /i 2 ). Doing the integrals, one 
obtains 

/L = ' J { 1 -^rr^)- m - (21) 

The result for n = is the same as in Q), while for n ^ it is new. To order 1/J, 
the vanishing points for higher pairs are the same as those for the lowest one. It must be 
remembered, however, that since we demanded x > 0, one must have A > 0. (Finding A < 
means that the oscillatory forbidden region has disappeared.) Thus the highest-field level 
crossing is successively eliminated as n increases, and (including zero and negative values) 
there are 2( J — n) fields in all where A vanishes. 
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FIGURES 

FIG. 1. Band-edge functions U±(m). For energy E a , mt and mi are turning points, the latter 
in the classically forbidden region. 
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